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Abstract 

We propose several T-odd asymmetries in the decay chains of the top 
squarks im — > and t — > bW~^ — > blu and Xk ~^ ^^^n ~^ l^l^Xi^ ^or 
I = e, fi, T. We calculate the asymmetries within the Minimal Supersymmetric 
Standard Model with complex parameters Mi , /i and At . We give the analytic 
formulae for the decay distributions. We present numerical results for the 
asymmetries and estimate the event rates necessary to observe them. The 
largest T-odd asymmetry can be as large as 40%. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) |2] with complex pa- 
rameters, there are new sources of CP violation in addition to the CKM phase of 
the Standard Model (SM). After redefining the fields these are the phase of the 
higgsino mass parameter fi, two of the phases of the gaugino masses Mi, i = 1, 2, 3 
(usually these are chosen to be the phases of Mi and M3), and the phases of the 
trilinear couplings Af, 0^^. The latter ones for the third generations / = T,t,b 
are rather unconstrained by the experimental upper bounds on the electric dipole 
moments of electron and neutron It is therefore especially interesting to search 
for observables which could be probed in forthcoming collider experiments in order 
to determine the phases (pA^tb- The influence of the phases 4>Artb been dis- 
cussed in the literature before. Some examples of studies discussing CP sensitive 
observables in SM processes or in processes which might occur in the SM with an 
extended Higgs sector are in jl] , where the influence of the phases (pAt t arises due to 
loop corrections. Other studies focus on the t b dependence in supersymmetric 
processes. There the dependence of (pA^ t b either CP-odd observables jH] or on 
CP-even observables jSl 13 El have been discussed. A CP sensitive asymmetry in the 
3-body decay ti hUrT"^ involving the transverse polarization of the r lepton has 
been proposed in [9j. 

In this paper we investigate whether the search for aplanarities in the decay 
chain of the top squarks ti^2 can give information on the CP phase ipAt or on other 
couplings of the MSSM Lagrangian. We consider the decay chain 

^ txl (1) 

with the subsequent decays of t and x^. We work in the approximation when both 
the top quark and the neutralino xl are produced on mass-shell. As the top quark 
does not form a bound state (because of its large mass), both the top quark and 
the neutralino Xfc decay with definite momentum and polarization. We consider two 
possibilities for the top quark decay: 

t bW+ and t bW+ blu{bcs), (2) 

and the following two-decay chains for xl'- 

xl^Ut. in ^12X1 and X°-C/r, in-^^tXi, (3) 

where the label of the leptons indicates their origin and where both if and are 
from the same lepton family. (Sometimes in literature if and are called the near 
and far lepton, see e.g. PUj.) We assume that the momenta of all ordinary particles 



2 



in ((H)-© can be measured or reconstructed, these are pt, Pb, Pi, Pi± and Pi^. The 
final state consists of two opposite signed leptons of the same family, a b quark 
and qq' jets (or /) from the t quark decay and missing energy. 

An useful tool for studying CP violation are triple product correlations (qi x q2 ■ 
qs) = (qiq2q3) ^H^l, where q^ can be any of the 3-momenta of the particles in the 
decay chain. Triple product correlations are an example of T-odd correlations that 
change sign under a flip of the 3-momenta q^ — q^ The time reversal operation 
T implies not only reverse of the 3-momenta and polarizations of the particles but 
also an interchange of the initial and flnal states. Because of the antiunitary nature 
of the time reversal operation, a non-zero value of a T-odd observable would imply 
T-violation if loop amplitudes are neglected. Any triple product correlation would 
be a direct evidence that T invariance is broken and as CPT invariance holds, CP 
conservation as well. As the triple product correlations in the processes (0)-® 
are a tree- level effect, they do not contain the suppression factor due to radiative 
corrections that is always present when such correlations are considered in processes 
with ordinary particles. 

In the top squark decays (H}-© no triple product correlations can arise solely 
from the decays of either t or xt- Triple products originate from the covariant 
products e{qiq2q3qi) written in the laboratory system. In order that e{qiq2qzqi) 7^ 
leads to a CP asymmetry at tree-level we need both a CP violating phase and at 
least a 3-body decay mediated by a particle that is not a scalar. The top quark 
decay modes © proceed in the SM and at tree-level no CP violating phases occur, 
thus no correlations of the type (pzP6Pt) can appear. The Xfc decays Q are 3-body 
decays mediated by the scalar lepton /„ and, as /„ does not transfer information 
about the spin of x\ to its decay products, again no triple products can be formed. 

Thus, the only correlations which occur are among the momenta of the decay 
products of both t and x\. These correlations reflect the spin properties of t and Xfc- 
In order to obtain analytic expressions for the distributions of the decay products 
we use the formalism of Kawasaki, Shirafuji and Tsai ^Hj- We work in the narrow 
width approximation for t and Xfc- 

As T-odd observables we consider up-down asymmetries, which are defined by 



^Jd^l sgn{0) dT/dn _ N[0 > 0] - N[0 < 0] 

/ dQ dV/dVL ~ N[0 > 0] + N[0 < 0] ' ^' 

where dV stands for the differential decay width and dVt involves the angles of 
integration. In Eq. O represents the triple product correlation on which we 
focus and N[0 > (<) 0] is the number of events for which O > (<) 0. According 
to the decay channels of the top quark we consider two cases: 
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1) If t — » bW, the possible triple products are 



(P6Ptlt2) and (pbKh)- 



(5) 



2) If a final leptonic (ul) or hadronic (cs) decay mode of W is measured, then 
possible triple product correlations are: 



In most of the asymmetries studied below 6-tagging will be necessary. In those 
asymmetries which involve the decay W —>■ cs also c-tagging will be necessary J3] . 

The decay xt ^ ^ l^l^ , leads to the same final state as decay (jHI) and 

gives also rise to the above triple product correlations. In this paper we will not 
consider triple product correlations in this decay, because due to the small Z^l^l~ 
vector coupling one can expect that the corresponding T-odd asymmetries are much 
smaller than those following from the decay 0- (In this context see also [131ll6| IT7j.) 

The paper is organized as follows: In the next section we give the relevant terms 
of the Lagrangian. In section El we present the results of our calculation in compact 
form using the formalism of Q^. Section |3] contains the formulae for various decay 
distributions. We propose several T-odd asymmetries in section El In section (HI 
we perform a numerical analysis of the T-odd asymmetries proposed. Finally, we 
summarize and conclude in section [71 



{VtViAt-i) ' iVbViAt^) and (p/,cl|l2 ) • 



(6) 



2 Lagrangian and couplings 



The terms of the Lagrangian necessary to calculate the T-odd asymmetries and the 
decay rates of im Xfc^ and Xfc In^f ~^ ^i^fXi in the presence of the CP phases 
are: 



9 Ha^k Pr + bnk Pl) xl In + h.c. , 



(7) 




9 1 {(^Ink Pr + blnk Pl) xl im + h.c. , 
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where Pl,r 
are defined 



|(1 =F 75), g is the SU(2)^ gauge coupling constant and the couplings 



as 



<k = {KiYAk + {K2rh' 



nk 



(9) 
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with 7l[^j being the scalar lepton mixing matrix and 



h'm = {hi,r = m,s, (10) 



and 





n=l n=l 

Here 7?.^„ is the mixing matrix of the top squarks and 

(12) 

\'^Rk/ \JRk/ 

with 

fik = -^{Nk2 + ltanewNki) , 

fkk = '^te^nOwK,, 

hik = ih'nkT = -YtN;, . (13) 

The unitary 4x4 neutrahno mixing matrix is defined in Appendix A, Eq. fj7ip . 
Yt = 'mt/iV^'nT'W sin (3) and Yi = mi/{\/2miY cos /?). The top squark mixing matrix 
7^* is given in Appendix B, the scalar lepton mixing matrix TZ^ can be found, for 
instance, in [H]. 



3 Formalism 

According to the formalism of the differential decay rate of (Q)-©, when spin- 
spin correlations are taken into account, is: 

dT = dT{t^ txl) dT{t ^ ...) dVixl -...), (14) 
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where the factors E^^/m^^T^^ and Et/rritTt stem from the used narrow width ap- 
proximation for if: and xty and T^^ are the total widths of the particles and m^^ 
and rrit are their masses. We have 

dr{t^-^txl) = ^\A\^d<l>i, (15) 



^77 



where 

is the differential decay rate of the top squark tm into a top quark with polarization 
4- vector and a neutralino with polarization 4- vector For the matrix 
element A we have: 

A = guipt){bl,PL + al,Pj,)v{p^,). (17) 
In evaluating |y4p we use the spin density matrices of t and xt'- 

P{pt) = MPt) P(-Pxj = -A(-PxJ ^^f^, (18) 

with 



HPt) =i>t + m^, A(p^J =i>^^ + m^^, (19) 

where "pt and p^^. are the momentum 4-vectors of the top quark and the neutralino 
X^. We have: 

„2 

i^r = y {(i«Lr+i&^„fcr)[(Px.pt)+'^x.'^t(^x.6)] 

^mfc) hxfc"^t - (Pxfe6)(^X.P0 + (Px.Pt)(^Xfe6)] 

+ 2 3m(aJ:,6^,)e(p^,ex.6p*)}, (20) 

where e^^^s _ rjj^g polarization 4-vectors and are determined through the 
decay processes of the top quark and the neutralino. dV{t ...) and c?r(x° — >• ...) 
are the differential decay rates of the unpolarized top and unpolarized neutralino. 

Next we shall consider the decays of and t. According to the chosen decay 
mode of the top quark, Eq. 0, we have to distinguish two cases. We consider them 
separately. 
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3.1 Decay rates for ifl^ 



For the width of the neutrahno decay into a lepton if and a scalar lepton /„ we 
write 

- m = ^:^Tr(5(-A(-p,J)i?)rf$,„ (21) 

Xfc 

with 

Here B is defined through the decay matrix element: 

%{p^,)B'^ = gv{p^M,Pn + J:kPLHlt), (23) 

p^f, and if are the momentum 4-vectors of the neutralino and the lepton, and E+ 
are the energies of Z„ and If. For the distribution of the decay products we obtain: 

dTixl - Ifl-) = Tp^iKkl' + \tk\')iPxJt) rf*x.- (24) 

For the polarization vector of the neutralino Xk^ determined through the Xk~ 
decay, we have: 



,a _l,a/3_ Px.PxA Tr (^(- A(-p^ J)757/3^) _ ^x. n- (9^\ 
<J Tr(5(-A(-p,J)i?) -""UPxjr^^' ^ ^ 



with 



Ql= ilfr-Si^llp^\^ «+= y' (26) 



Respecting the condition (^^fc^'xfe) ~ 0' vector is orthogonal to the momentum 
4- vector of xt- This is the only orthogonal 4- vector composed of the available 
momenta p^^. and If. As it can be seen from and is in the x^-decay 

plane. Further we shall assume that In is produced on mass-shell, pj = rrij, then 
{PxJf ) = ("^Xfe ~ ''^[)/2) where we neglect the lepton mass in the kinematics, i.e. 
mi = 0. The prefactor a+ determines the sensitivity to the polarization of xt- 

Note that the polarization vector of xl, Eq. (|^. does not change if we take the 
subsequent decay /~ — > Xi^2 i^ito account. Note further, that the polarization vector 
of the C-conjugated decay Xk ~^ ^i ^n changes sign compared to ^^j. in Eq. 
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3.2 Decay rate for t bW+ 

When the top quark decays according to t — bW^ we have: 

1 



dT{t ^ bW+) = Tr(aA(Pt)a) d<!>1 , (27) 



with 



where Cb is defined by the decay matrix element as follows: 

C^u^ipt) = uipbh^PMPt) e'-^Pw), (29) 

where pb, Pt and pw are the momentum 4- vectors of the bottom quark, the top 
quark and the W boson. Then for the distribution of the decay products we obtain: 

8Et rriy^r 

We denote by ^b the polarization 4- vector of the top quark, determined by the decay 
t — > bW. Its expression is given by the formula: 



.a _ I „at3 _ PtPt \ Ti{CbK{pt)i^-ipCb) , . 

^^"'^ Tr(aA(j>t)a) ' ^ ^ 



From ()29|) and (jM^ we obtain the polarization vector: 

(32) 



with 



"^= mf + 2m^ ' (l'.P.) = ^^, ^3. = (^P.--I^P. h (33) 



where in the kinematics we have set rrib = 0. Here Qb is the 4-vector orthogonal to 
pf and ab determines the sensitivity to the polarization of the top quark. 
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3.3 Decay rates for t hlu and for t bcs 

We consider here the decay t — >■ blu. For the decay t — >■ bcs one has to make the 
replacements p^, — > Ps,Pi Pc in the equations below. For the inclusion of QCD 
corrections to the decay of a polarized top quark we refer to |TH]. When the top 
quark decays according to t hlv we have: 

dV{t ^ blv) = ^-1^ TT{CMPt)Ci) d^l , (34) 

with 

where Ci is defined through the decay martix element as follows: 

C^lUM) = {A^)MPuhaPLv{pi) ^""^ - P^^^/"^^ u{p,hpPMPt), (36) 

with 

Dw = (Pw - ^w) + irnwTw, Pw = Pt - Ph ■ (37) 
Then we obtain: 

rfr(t ^ Uv) = f (pm) {m\ - 2{pm)) d^l. (38) 
From ()36|) . for the polarization vector of the top, that we denote by ^i, we have: 

(„af3 P?Pt\ Tr(QA(pt)757/3C/) , m 



with 



Qr=br-^p?V «^ = -i, (40) 



mt 



Qf is orthogonal to and lays in the top quark decay plane (in the rest frame of 
the top quark). In general, with the available vectors in the decay, pt, Pb and pi, one 
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can form three independent combinations orthogonal to : two in the decay plane, 
and Qf, and one transverse to it, e{aptPbPi)- As CP invariance holds in the top 
quark decay, there is no contribution to the transverse component, and because of 
the V-A structure of the interaction, there is no contribution to either. 

Inserting and (I27j) or into jH)) we obtain dT in terms of the 

polarization vectors: 

^r''-' = / V |ApTr(^(-A(-p,J)i?)Tr(a,,A(p,)a,0 

X ^/(l«iir + l&iiP)K-Orf<J^"\ (41) 

where we have used the narrow width approximation for the scalar lepton propaga- 
tor, mi and is the mass and the total decay width of /, |Ap is given by pOjl and 
(l^b,i (jgnotes the phase space for the two different decay modes of the top quark: 

^(^b,i ^ . ^^b,i . ^^^^ . (^42) 

In order to obtain the angular distributions of the ordinary particles in (|T|) - we 
have to use the explicit expressions for and C,t, and carry the integration over 
the phase space of the supersymmetric particles. 



4 Decay distributions 

In this section we derive the analytical expressions for the decay distributions of 
im iXk ~^ hW^ 'X.i^t^2 and im — ^ tXfc ~^ bluxllfl^. We consider separately the 
two decays (j2I) of the top quark. 

4.1 Decay distribution for tXk bW^ Xi^t^f 

We choose in the direction of the Z-axis and and pf, determine the YZ-plane: 

Ptr. = Pt = IpI (0,0,1), Px, = IpI (0,0,-1), 

(43) 

where we have used the brief notation Cb = cos 9b, S(j,^ = sin0+, etc. The ranges of 
the angles are < 6b, 6 +,6- < vr; < 0_ < 27i. Then we can carry out part of 
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the phase space integration. Using (jinj) . (ji^ is given by 



, ^ IpI (^i - m^w) i^l, - ^f) dcb dc+ d(j)+ 



where 



1 rrir — ml 

d^i = ^7^2 r.2n a^\2 ' (45) 



and 



8(2vr)2E?(l-A-c,v)^ 



. _ IpI - IpI r - Ip^I i„i - ^'^'(^L>^oO 



2 I m2 _ tm * Xfe 



2mf 
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^x. = VlPl' + <= 2m- ^' ^«- = (PrP.-), (46) 



where A(x, y, z) = + y'^ + — 2{xy + xz + yz). Then from (HH), using the exphcit 
expressions for the polarization vectors ()25|) and ()32p and the decay distributions 
and (jnni), for the angular distributions of the 6-quark and the leptons li and 
Z^, we obtain: 



dch dVL> dQ_ 



^{{\<k\ +Wrnk\) 



iPtVb) {PxJi 



(|a: 



t |2 



mfc I J 



m 



7T (QbVxk) + «+ 



Xfe 



(Q+Pt) 



X 



(PtPb) {PxJl 



(PtPb) {PxJi 

liPxkQb) (Q+Pt) - (PxkPt) (Q+Qb)] 



+2 ^m{a^J}^^) ab a+ 



rat 'm 



Xk 



{PtPb) {PxJl) 



(l^PfcPt) 



(47) 
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where 



m^-i I S Air ^ 



a^Y (2m^ + m,^)K - ^wfi^, ' ^1? (Wnk? + 1^1') |p| 
27r m}^ rut Ttm^^ T^^ E^^ 



ml -ml, m\^ - m] - '^xi 

tLih = _ „ r, -C/4- = ^r^: 7-. ^ r, -C/- ~ 



2E,(1-Ac.)' " 2E,,(l + /5,,c+)' - 2Ej{l - Pjcj^.y 
'itpbPt) = E+Eb\p\sbS+c^^. (48) 



4.2 Decay distribution of txl blh'Xiltl2 

The angular distribution of the final 6-quark and leptons /, If and is obtained 
from the previous results if we fix the coordinate system so that pt and p; determine 
the YZ-plane: 

Pt = IpI (0, 0, 1), p^, = IpI (0, 0, -1), pi = El (0, Si, q), 

Pb = Eb (sbC0,, SbS^^, Cb), (49) 

where the ranges of the angles are < 6';,, 6'; < vr; < 06 < 2tt. Then the dependence 
on the 6-quark momentum is only in the phase space. We obtain: 



IpI {ml - m1^){ml^ - mf) 
2mi^8^i2nyE?El 

dci dVtb dVL^ dsw ,^ 

(1 - Ac6)2(l + /3^,c+)2[E,(l - Aq) - E,{1 - cuW ' ^ ' 

(50) 

where sw = Pw- obtain the angular distribution by a replacement of the phase 
space d^^ —>■ and the following replacements in the curly brackets of (jTTj) : 
ab ai, Qb ^ Qi and pb ^ Pi- The angular decay rate distribution of Zj*", Z^, I and 
b is: 

Nr- 



dci dQb dn+ dn_ ^ ' (1 - (3j )2 



I 

{.PtPi){ml - 2{ptpi)) 



(1 - Ac,)2(l + p^^c+f[Et{l - Aq) - Eb{\ - Cbi)f 

X l(i«Lr+i&Lr) 



(p,,p.)+a,a.^^(Q,Q,) 
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- a: 



t |2 
mk I 



rmfcl J 



[PtPi) [PxJi 



- 2 3fJe(a^fc6^fc)m^^mt 



X 



1 - 



(PtPl) {PxJl 



[{PxkQi) (Q+Pt) - (PxkPt) iQ+Qi)] 



+2 $5m(a46„fc) ai a+ 



rrit m 



Xk 



(PtPl) (PxJl 



(51) 



where 



N, = 



El 



m 



w) \P\ 



Xk 



(|a: 



nk I 



\"nk\ ) 



nit ^tm^k Txfe "^vy 



Xk 



m 



w 



2[Et{l-f3tCi)-E,{l 
E+Ei\p\sis+c^^. 



Cbl)] ' 



Cbl 



[Pb ■ Pi) 



(52) 

(53) 
(54) 



The distribution of and I2 from the C-conjugate decay xt ^ h^t ^ h ^2X1 is 
obtained from ()47j) and 1)511) by the replacements if li , I2 it — —«_)_. 



As can be seen from the angular distributions, Eqs. and ()51|) . the prefactor 
of the triple product correlations (last term in Eqs. (jlTj) and (jSH)) is ^''^(amfc^mfc) 
and consequently the T-odd asymmetries (to be defined in the next section) are 
propotional to this prefactor. Therefore, in order to study the dependence of the T- 
odd asymmetries on the MSSM parameters, it is useful to give the explicit expression 
for ^mi^a^^jj^l^^) for m = 1 using Eqs. (fTT|) - (fT^ : 



53m(a*;;6*,) = - cos^ ^,-Fi53m(/i/iV*J - sin^ tan ^^53m(iV*,iV,*J 

+ cos sin % (^^tan^H^53m(/i/iV*,e^<^^) + Y^^m{Nl,Nl,e~'^-^)^ . 

(55) 



We can see from that if CP violation is solely due to (pAt 7^ 0, the T-odd 
asymmetries are proportional to sin 2di sin (pi, which can be naturally large because 
of the large top squark mixing (see Eqs. ()72|) and (|73j) l. Moreover, one can see from 
()55|) that the term oc sin26'j can be sizable also in a higgsino-like scenario (|yu| < M2) 
because of the large top Yukawa coupling. 
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5 T-odd asymmetries 



We shall distinguish three classes of asymmetries according to the lepton momentum 
(stemming from the decay chain x° l^lf ^t^2^Xi) involved in the triple product: 
(i) when the momentum vector of lepton if from the decay ^n^t enters; (ii) 

when the lepton momentum vector of lepton /J from the decay — > /Jxi enters; 
and (iii) when both momentum vectors of if and /J from the decay xt ~^ ^JXi 
enter. 

• The first class involves the asymmetries: 



N 


(Pbpaf ) > 


-N 


'{pbPtlf) < 


N 


'{PbPtlf) > 


+ N 


'{PbPtlf) < 



(56) 



N 


\piPtlf) > 


-N 


'(piPtlf) < 


N 


'(piPtlf) > 


+ N 


\piPtlf) < 



(57) 



N 


(p^p,l±) > 


- 


'{piPblf) < 


N 


(piPblf) > 


+ N 


(p^paf ) < 0^ 



(58) 



where p; is the lepton momentum in the decay t — > blv. 

• In the second class of the asymmetries h is replaced by I2: 



N 


(pbPtl^) > 


-N 


(PbPAt) < 


N 


(P6P*1^) > 


+ N 


\pbPtit) < 



(59) 



N 


\piPtlt) > 


-N 


'{piPtlf) < 


N 


\piPtlt) > 


+ N 


\piPtlt) < 



(60) 





(piPblt) > 


-N 


(piPblt) < 


N 


(P;P&1^) > 


+ N 


(piPblt) < 



(61) 



• The third class of asymmetries is: 



N 


\pblf\^) > 


-N 


\pbifm < 


N 


'(Pblfl^) > 


+ N 


\pblflf) < 



(62) 
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At 



N 


baf ID > o; 






N 




+ N 





(63) 



Since the polarization vectors of xl for the two C-conjugate decay modes of ^ differ 
only by a sign (see Eq. (1221)) ^^e value of the asymmetries with upper indices + and 
— are related by: 



A1 



-A: 



,5) 



and 



-A 



1,2,3). 



(64) 



In order to measure all of the listed asymmetries it is necessary to distinguish the 
lepton if, originating from the decay Xfc ^n^t^ the lepton from the subse- 
quent decay Xi^2- This can be accomplished by measuring the energies of the 
leptons and making use of their different energy distributions, when the masses of 
the particles involved are known, if and l^ can be be distinguished if their measured 
energies do not lie in the overlapping region of their energy distributions. 



We define the fourth class of asymmetries as follows: 
N [(p,pa+) >0]-N [(p,pil+) < 0] 



A, 



A, 



A, 



iV[(p,p,l+) >0]+iV[(p5Pil+) <0]' 

iV[(p^p,I+) >0]-iV[(p,p,l+) <0] 
iV[(p,pil+) >0] + iV[(p,pil+) <0]' 

N [(pipa+) >0]-N [(pipa+) < 0] 



(65) 



(66) 



(67) 



N[{pip,l+)>0]+N[{pip,l+)<0] 

where l"*" stands for the momentum 3- vector of either If or Evidently we have 

Af + Aj'^ 



A, 



1,2,3. 



(6^ 



A measurement of Ai, i = 1,2,3, does not require to distinguish between the leptons 
if and 12^, it requires only a measurement of their charges. Analogous asymmetries 
can be defined for /~ as well. 

It should also be noted that the asymmetries given in Eqs. (fH m) - (j3H|) do not 
depend on the mass of In- The asymmetries above are written down for the leptonic 
decay lu. For the hadronic decay — > cs the analogous asymmetries are 

obtained by replacing p^ — > pc. 
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6 Numerical results 



All proposed T-odd asymmetries depend on Qm{al^i.bl^i^) , Eq. (j^ . and measure 
therefore the same combination of CP phases in the MSSM, but they have differ- 
ent magnitude. In this section we present numerical results for the asymmetries 
Af,A'^ {i = 1,2,3), Eqs. lEB-^, ^J,^5, Eqs. ^ and ^ and Ai {i = 1,2,3), 
Eqs. fl(i5j) - (jF)7jl . In order not to vary too many parameters we fix m^^ = 400 GeV, 
771^2 = 800 GeV and tan/? = 10^. We take rrit = 178 GeV and we also use the 
GUT relation \Mi\ = 5/3 tan^ Bp^Ms. We take m^-^ = 130 GeV, mj^ = 300 GeV, 
where we assume li ~ Ir for I = e, n, which is suggested in mSugra models. In the 
scalar tau sector we take into account scalar tau mixing choosing A^ = 500 GeV. In 
our numerical study we take |y4t|,0^j, M2, (pMi, \^J'\,4>^J. as input parameters. Note 
that for a given set of input parameters we obtain two solutions for (Mg, M^) cor- 
resonding to the cases M^^ > Mf^^ and M^^ < M^^ in Eqs. ^ and (HH) which 
we will treat separately. In the plots we impose the phenomenological constraints: 
m-± > 103 GeV,mj^o > 50 GeV and Xi is the lightest supersymmetric particle 
(LSP). 

In Fig. ^we plot the various asymmetries ()56 |) -()67 |) for the decay ti a 
function of (pAt for the case Mq < Mjy. The MSSM parameters are M2 = 250 GeV, 

= 200 GeV, \At\ = 1200 GeV and 0Mi = 4>f, = 0. As can be seen in Fig. ^ the 
absolut value of the asymmetry A2 (dashed line) is much larger than the absolut 
value of Ai (solid line), which can be attributed to the sensitivity factor of the top 
quark polarization for the two asymmetries. For Af this factor is \ab\ — 0.4 (for 
mt = 178 GeV), Eq. whereas for it is |az| = 1, Eq. (gU)) (see also [iS] 

where QCD corrections are included). This difference can also be seen in Fig. ^ 
by comparing A^ , Eq. (jHH), with A^ , Eq. (jUHj) . In Fig. [TJi the asymmetries Ai {i = 
1,2,3), Eqs. (jniI)-(|nZj)! are displayed for which we have to distinguish the leptons 
in the decay chain X2 ~^ ^f^t ~^ Xi^J^f only by there charge. It is interesting to 
note that the asymmetry A2 can be as large as 24%. In Fig. |21 we plot the same 
asymmetries as in Fig. Q but now for the case Mq > M^. As can be seen also for 
this case the largest asymmetry is A2, Eq. ^7\ . which is however somewhat reduced 
compared to the case Mq < Mfj. 

In Fig. 121 we plot the contours of the asymmetry A2 for the decay ti tx2 in 
the |/i| — M2 plane where we have taken (pAt = f and the other parameters as in 
the previous figures. One sees in Figs. 12^-1211 that the asymmetry is largest for large 
gaugino-higgsino mixing {\fi\ ~ M2). Figs. 12^' and b correspond to the case where 
X2 decays into I = e,fi, wheras Figs. |2t and d correspond to the case X2 ~^ tit"*". 
The asymmetries in Figs. |2t' and b are larger than those in Figs. |2t and d, because 

^The prefactor Srn(a^j,6^jj,) of the asymmetries does not depend very much on the value of 
tan/3 if \At\ > |/i|/tan/3. 
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4>At I TT 



Figure 1: T-odd asymmetries (a) v4+, i = 1,2,3, Eq. (b) A'r ^ i = 1,2,3, 

Eq. dSnil-dnil), (c) Af, i = 4,5, Eqs. (jSl and and (d) Ai, i = 1,2,3, Eqs. 

for ii as a function of 0^^. In (a), (b), (d) the solid (dashed, dotted) hues 

correspond to the indices i = 1 (2,3), in (c) the sohd (dashed) hne corresponds to 
i = 4 (5). The MSSM parameters are chosen as |A| = 1200 GeV, Ms = 250 GeV, 
= 200 GeV, tan (3 = 10, = (j)f, = 0, mi^ = 400 GeV, = 800 GeV, 
Mr) < M^, for l = e,iJ,. 



of the effect of scalar tau mixing which leads to |a+| < 1 (see Eq. ()26|) ). while for 
I = e,fi, Figs. ^ and b, we have \a^\ = 1. Moreover, in Figs. and b there is a 
sign change of the asymmetries because of a sign change of the prefactor 53m(a^2^5_2)- 
This sign change does not appear in Figs. Et and d, because it is compensated by 
a simultaneous sign change of a+, which occurs due to a level crossing of the states 
X2 and Xs- Figs. Et and d there is a different sign change of the asymmetries in 
the lower right part of the — M2 plane because a+, Eq. (|^. changes sign there. 

In Fig. m we show the contours of the asymmetry A2 for the decay ti 
bW+lilt {I = e,/i) in the 0^* - 0Afi plane taking At = 1200 GeV, M2 = 250 GeV, 

= 200 GeV and 0^ = for the two cases Mq < Mq (Fig. Hi) and Mq > Mfj 
(Fig. IUd). In Fig. the largest value of about 39% (—39%) for the asymmetry 
A2 is obtained for (pAt = l-47r (0At = O.Gvr). For (pAt = (tt) and = 0.57r 
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Figure 2: T-odd asymmetries (a) Af, i = 1,2,3, Eq. (b) A'r ^ i = 1,2,3, 

Eq. dSnil-dnil), (c) Af, i = 4,5, Eqs. (jSl and and (d) Ai, i = 1,2,3, Eqs. 

for ii as a function of 0^^. In (a), (b), (d) the solid (dashed, dotted) hues 

correspond to the indices i = 1 (2,3), in (c) the sohd (dashed) hne corresponds to 
i = 4 (5). The MSSM parameters are chosen as |A| = 1200 GeV, Ms = 250 GeV, 
= 200 GeV, tan (3 = 10, = (j)f, = 0, mi^ = 400 GeV, = 800 GeV, 
Mq > M^, for l = e,iJ,. 

the asymmetry A2 is about —3.5 (17.3)%. In Fig. one can see that the largest 
value of the asymmetry A2 is obtained if (pMi 7^ 0, tt and the asymmetry varies from 
about 25% for (pAt ~ I-Stt, 0Mi = to about 34% for (pAt ~ I-Stt, ~ lAn. For 
(pAt = (tt) and 0Mi = 0.57r the asymmetry A2 is about —10.6 (3.2)%. 

In Fig.Elwe show the asymmetry A2, Eq. (jSZI), as a function of (pAt ioi the decay 
of the heavier top squark. The MSSM parameters chosen are \At\ = 1200 GeV, 
M2 = = 200 GeV and (pMi = (pii = 0. Fig. displays the asymmetry for 
the decay ^2 — ^ ^Xs and Fig. shows the asymmetry for ^2 — ^ ^X4- Figs. and b 
demonstrate that also for the decay of ^2 into the heavier neutralinos the asymmetry 
A2 can be quite large. In Fig. the two cases Mq < Mjj and Mq > Mq- give nearly 
the same curve for A2, because the value of the decay width is nearly the same for 
the two cases. Note that an observation of a T-odd asymmetry would lead to a two 
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Figure 3: Contours of the T-odd asymmetry A2 in % for ti ^ bW'^li if for 

I = e,fj., (a) Mq < and (b) Mg > M^j and for / = r (c) Mq < and (d) 
Mq > Mfj. The MSSM parameters are chosen as \At\ = 1200 GeV, (pAt = O.Stt, 

400 GeV, mr = 800 GeV, \AJ = 500 GeV, 



tanp = 10, 0Mi = 0^ = 0, m^-j = 
4>Ar = 0, nif-^ = 130 GeV and m^^ = 300 GeV. The hght gray region is excluded 
because there m^± < 103 GeV and/or m^o < mj_^. In the dark gray area the 

two-body decay ii — > is kinematically forbidden. 



fold ambiguity in the extraction of the CP phases, which can be seen in Figs. I1|2|4I 
and El 

Next we give a theoretical estimate of the number of top squarks ti necessary to 
observe the T-odd asymmetries (f^ - (pT7j) in the decay ii — > X2^y where subsequently 
the neutralino X2 decays into lil (/ = e,fi,T). This number can be estimated by 



Nr > 



a 



(69) 
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4>At/T^ 



Figure 4: Contours of the CP asymmetry A2 in % for ti —>■ ^ bW^li l^. The 
MSSM parameters are chosen as \At\ = 1200 GeV, M2 = 250 GeV, = 200 GeV, 
- w, - 10, rrii^ = 400 GeV, m^^ = 800 GeV, for I = e, fx, (a 



0, tan/3 
and (b) Mq> M^. 

(a) 





Figure 5: T-odd asymmetry A2, Eq. (jSTj), for the processes (a) t2 tXs ~^ bW'^li If 
and (b) ^2 — ~^ bW~^lili (/ = e, yu) as a function of (pAf The MSSM parameters 
are chosen as |A| = 1200 GeV, M2 = = 200 GeV, tan/? = 10, (j)Mi = 0^, = 0, 
m^-^ = 400 GeV, m^^ = 800 GeV for Mq < (sohd hue) and Mq > (dashed 
line) . 



where a denotes the number of standard deviations and At stands for any of the 
above asymmetries. The value of the branching ratios of the W boson is given 
by B{W ^ f) = (32,68,32)% (/ = Y.ivl.Y.qqq\cs) HH] corresponding to the 
asymmetry for which the estimate is made. For instance, for the asymmetry , 
which is based on the triple product (pbpjf ), we take the value for B{W Y.q Qq')- 
The estimate is taken for two scenarios which we define in Table ^ The scenarios 
chosen imply that J2i B{x2 ~^ h^) = 1 (^ = e, /i, r) and that the energy distributions 
of the two final leptons (from the decay xl ~^ ^1^1 from the decay — >• X1I2) do 
not overlap. This means that in all decays it is possible to decide from which decay 
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Table 1: The two scenarios used for the estimate of the necessary event rates. 



scenario 1 


scenario 2 


= 129 GeV 
M2 = 500 GeV 
= 150 GeV 

tan (3 = 3 

\At\ = 1200 GeV 

= f 
Mq<M^ 


mj^ = 115 GeV 
M2 = 200 GeV 
1^1 = 300 GeV 
tan/5 = 6 

4>Mi = 0, 0^ = TT 

\At\ = 1200 GeV 



Table 2: The values of the T-odd asymmetries defined in Eqs. (fSU jl - ljUTI) and the 
number A^^^ of top squarks required to measure these asymmetries with a 3a evidence 
in the two considered scenarios (see Table HI). 





scenario 1 


scenario 2 


At 


value [%] 


Ni, ■ 10-3 


value [%] 


Ni, ■ 10-3 


At 


-11.5 


4.5 


15.9 


24.8 


At 


28.3 


1.6 


-39.0 


8.7 


At 


13.8 


6.8 


-16.3 


50.4 




-1.3 


355.9 


5.1 


242.2 




3.2 


124.9 


-12.6 


84.3 


A's^ 


1.6 


499.2 


-5.3 


476.8 


At 


-4.7 


18.5 


6.1 


115.0 


At 


11.4 


9.8 


-14.9 


60.2 


Ai 


-6.4 


14.8 


10.5 


57.1 


A2 


15.8 


5.2 


-25.8 


20.1 


A3 


7.7 


21.6 


-10.8 


114.8 



the two leptons originate. For simplicity we will assume that the T-odd asymmetries 
(|3H|) - (jU7|) are equal for the three flavors in the subsequent decay X2 ~^ h^- This 
means that we neglect scalar tau mixing (|a+| = 1 for / = e,ii,T) and in addition 
we take mj_^ = mg^ = m^^ = m^^. For the calculation of the branching ratios 
of the ti we use the formulae given in [7j. For scenario 1 (scenario 2) we obtain 
B(ti X2t) = 22% (2.1%) where we assume that the bosonic decays ti biW^ 
and ti — s> biH^ are kinematically not accessible. 

In Table 121 we display the values of the asymmetries, Eqs. (j^ - (jU7jl . and the 
numbers Nf^ needed for a 3a evidence of these asymmetries. From Table |21 it can 
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be seen that in order to have a 3a evidence for some of the T-odd asymmetries in 
scenario 1 about O(IO^) produced ti's are necessary. For scenario 2 0(10"^) produced 
ti's are necessary for a 3a evidence of some of the T-odd asymmetries. Assuming 
that O(IO^) ti's can be produced at the LHC and O(IO^) ti's at a future hnear 
coUider, there are good prospects to measure some of the asymmetries. It is however 
clear that detailed Monte Carlo studies taking into account background and detector 
simulation are necessary to predict the expected accuracy. This is, however, beyond 
the scope of the present paper. 

7 Summary and conclusion 

We have proposed a set of T-odd asymmetries in the decay im — ^ txt with the 
subsequent decays t hW^ blu and xl 11^ I'^PxiJ for ^ = e,fi,T. The 
asymmetries are based on triple product correlations involving the polarizations of 
the top quark and the xt ^"^^ arise already at tree-level. All the proposed T-odd 
asymmetries probe CP violation in the t — tm — Xk couplings and are proportional to 
the product of left- and right-couplings. Since top squark mixing is naturally large 
due to the large top Yukawa coupling these asymmetries may be large and will allow 
to determine the CP violating phase (pAt, which is not easily accessible otherwise. 

In a numerical study of the T-odd asymmetries we have found that the asym- 
metry Af, which is based on the triple product {pipdf), is the largest one and its 
magnitude can go up to 40%, while the others are smaller. We have also found that 
the asymmetry A2, Eq. ()fi5|l . based on (p^ptl^), where can be any of the final 
leptons if and /J, distinguished only by their charges, is < 26%. Moreover, we 
have made a theoretical estimate of the number of ti necessary to observe the T- 
odd asymmetries for two scenarios. Depending on the MSSM parameters, we have 
found that a ii production rate of O(IO^) may be sufficient to observe some of the 
proposed T-odd asymmetries, which could be possible at the LHC or at a future 
linear collider. 
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Appendix 



A. Neutralino Masses and Mixing 

At tree-level the neutralino mass matrix in the weak basis {B, W^, H^, H^) is given 
as 



M 



N 



( iMile*"^^! 

M2 

-mzswcp mzcwcp 

\ rnzSwS/B -mzCwSp 



-mzswcp mzSwSfi \ 
mzCwCfs -mzcwsp 






— l/ile**^" 



— |/i|e*'^'' 



(70) 



where 0Mi is the phase of Mi, and cw and sw are cos 9w and sinOw, respectively. 
This symmetric complex mass matrix is diagonalized by the unitary 4x4 matrix 
N: 



N*MnN'^ = diag(m^o, . . . , m^o) 



< m^o < . . . < m^o . (71) 



B. Masses and mixing in squark sector 

The left-right mixing of the top squarks is described by a hermitian 2x2 mass 
matrix which in the basis (tL,^ij) reads 



tLL 



R) 



2 I 



tRR 



(72) 



where 



Ml 



M? 

tRR 



tRL 



M| + (- - - sin^ Qw) cos 2/5 m| + ml 



2 3 



2 

Mj^ -|- - sin Qw cos 2/3 + , 



(M|^ J* = mi(A-/x* cot /3), 



arg[At -/i* cot 



(73) 

(74) 
(75) 
(76) 
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where tan/5 = V2/V1 with Vi{v2) being the vacuum expectation value of the Higgs 
field Hf{H2), nit is the mass of the top quark and Qw is the weak mixing angle, /i is 
the Higgs-higgsino mass parameter and Mq, Mjy, At are the soft SUSY-breaking pa- 

rameters of the top squark system. The mass eigenstates U are (^1,^2) = (tL, t^jTZ^ 
with 



( e**^' cos Qt sin Qi 



\ — sin dt e '"^f cos Qi 



1 



(77) 



with 

-|M^ I M? -mj 

cos 9i = sin 9i = .m 

J\M? |2 + (m? -M? )2 ,/|M| |2 + (m? - M| )2 

The mass eigenvalues are 
1 



TTlr 



Note here that for \At\ ^ cot/9 we have 
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